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Abstract 

This paper is about the rate of convergence of the Markov chain X n +i = ylX„+B„ 
(mod p), where A is an integer matrix with nonzero eigenvalues and {B„} ?1 is a 
sequence of independent and identically distributed integer vectors, with support not 
parallel to a proper subspace of Q fc invariant under A. If |Aj| ^ 1 for all eigenvalues 
Xi of A, then n — O ((lnp) 2 ) steps are sufficient and n = 0(\np) steps are necessary 
to have X n sampling from a nearly uniform distribution. Conversely, if A has the 
eigenvalues Xi that are roots of positive integer numbers, |Ai| = 1 and |A; > 1 for all 
i y£ 1, then O (p 2 ) steps are necessary and sufficient. 
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1 Introduction 



In this paper we generalize some results obtained in the paper [2] about Markov chains on 
Z fc of the form 

X n+ i = AX n + B n (mod p), (1) 

where Xo = xo G Z k , A € GLfc(Q) n M^.(Z), p is an integer and {B„,} n is a sequence of 
independent and identically distributed integer vectors. 

If k = 1 and B n is a fixed integer b, for particular values of p this recursion is used to 
produce pseudorandom numbers on computers (see, for example, Knuth's book [9]). 

In the paper [3], the constant term b is chosen with a fixed probability at each step 
and the authors study the following Markov chain: 

X n+ x = aX n + B n (mod p), 

where a is a positive integer. This randomness is introduced in order to produce uniformly 
distributed random numbers on the set {0, l,...,p — 1}. In the cited paper, it is shown 
that, for a = 2, n = 0(lnp In hip) steps are sufficient to sample X n from an almost 
uniform distribution. On the other hand, if a = 1 then n = O (p 2 ) steps are necessary 
and sufficient to achieve randomness. A further generalization is described in [6|, where 
the integer a is also allowed to vary. 

The extension of the previous results to the higher-dimensional case is due to Asci 
[2] and next to Hildebrand and McCollum [8], with the study of some particular cases of 
the recursion ([T]). In [2], the distribution of B n is the most general (the support of the 
distribution cannot be parallel to any proper subspace of Q fc invariant under A), but the 
matrix A has only integer eigenvalues. In [8], A is arbitrary, but only a specific distribution 
for B n is considered. 

The general case is studied in this paper, with the condition |[B n || 00 £ L 2 and some 
further conditions on p. We find two different types of behaviour for the sequence ([T|), 
depending on the size of the complex eigenvalues of A. If |Aj| ^ 1 for all eigenvalues Aj, 
then n = O ((lnp) 2 ) steps are sufficient and n = 0(lnp) steps are necessary to reach the 
uniform distribution (theorems 3.1 and 3.13). In particular, for a matrix A with eigenvalues 
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Aj that are roots of positive integers and |Aj| > 1, we can show that n = O (In pin In p) 
steps are sufficient (theorem 3.7). On the other hand, if the eigenvalue A, are roots of 
positive integers, |Ai| = 1 and |Aj| > 1 for all i ^ 1, then O (p 2 ) steps are necessary and 
sufficient (theorems 3.12 and 3.14). These theorems agree with the one-dimensional case 
studied in [3] and with the results in [2] and [8]. 

In Section 2, we recall some general results about random walk on groups and the 
preliminary lemmas proved in [2]. The main results of our work can be found in Section 
3. 

2 Preliminary results 

Consider the sequence (P) and observe that we can suppose X n G Z^. 

Set P n (x) = p(X n = x), V x G Z*, and /i(x) = p(B n = x), V x G Z k , V n G N; 
moreover, denote by U the uniform distribution on Zp. Define: 

V = {x G Z fc : x = h — k, where h, k G supp ^}. 

Indicate by d, where d < k, the degree of the minimum polynomial of A. By definition: 

d d 

Y[(A - Xil) = A - XJ) = G Mfc(Z), Aj G {Ai, X d }, V i = d+ 1, fc, 

8=1 1=1 

where Ai, Ad, Afc are the eigenvalues of A. Finally, set: 

V d-i = |^m x :x€ y ,m = 0,l, d - 1}. 

In order to show that the distribution P n tends to the uniform distribution U, as 
n — > +oo, use the Fourier analysis (see Diaconis' monograph jl] and [5]). Define the 
variation distance between P n and U in the following way: 

11^-011 = ^ E \Pn(a)-U(a)\. 

It is possible to prove that 

||P n - U\\ = \ sup \E Pn (f) ~ E v (f)\ = max \P n {A) - U(A)\, 
2 f eF AcZ£ 
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where F = {f — ► C :||/|| < 1}. 

Henceforth, our purpose will be to find a bound for \\P n — U\\ in terms of n and p. 
Observe that, if we indicate with {Y n } n the sequence defined by ([T]) and the condition 
Xo = 0, we have X n = ip n (Y n ), where the one to one function ip n : — > is defined 
by </? n (x) = A n XQ + x. Moreover: 

\\Pn-U\\ = \\(P n o <p n ) - U\\, 

then we can consider Xo = 0. 

Let / : Z fc — > C; define the Fourier transform / : R fc — ► C by: 



/(a) = Y, exp(^(h,a)) /(h). 



hez k 

We have the following four results, whose proofs are similar to the proofs of the lemmas 
2.5, 3.1, 3.3 and 3.4 in [2]: the only difference is that a ranges in R fc instead of in Z^. 
Lemma 2.1 is also proved in [1], in a more general case. 

Lemma 2.1 (Upper bound lemma). 

\\Pn-Uf<\ £ |P„(a)| 2 . 

a£Z£-{0} 

Lemma 2.2. Suppose gcd(det(^4),p) = 1, Xo = 0, a G R fc ; then: 

n-1 

1) Pn(a) = [] C^'a) . 

2) l p ™(«)| 2 = n(| E MhMi)cos(^<h-i,^a))j 

< I] (l - 2/i (u) /x (v) + 2/i (u) /x (v) cos <u - v, *^'a>^ ) , 

Vu, v G supp [i. 

Lemma 2.3. Let aeZp- {0}; then: 



\Pn~U\\>\ 



PJa. 



Lemma 2.4. Suppose that the support of \i is not parallel to a proper subspace of 
Q fc invariant under A. Then, there exists a basis {yi, •••,yfc} C V d ~ x of Q fc . Furthermore, 
for all p £ N such that gcd(det(yi, yk),p) = L for all a G R fc — (pZ) k , there exists 
i € {1, k} such that (yi, a) 7^ mod p. In particular, if the support of fi is not parallel to 
a proper subspace of Q fc , we have yi, ...,y& € V, (yi, a) 7^ mod p, for some ? G {1, k}. 

Henceforth, we will indicate by B the matrix (yi--.y&), where the vectors yi, y/c are 
defined by Lemma 2.4. 

Lemma 2.5. Vej'G N, we have: 



e-l 



i=l s=0 \ i=l 



^ e = H(U - A,/) + J; \ s+ iA e ~ s HC A - Xil) . 



*A j fie a - Xii) = j2 e n A - _m+i n ? a - a - j ) - v). 

i=l h=e+l ki,...,kd—h+l' m=h n=h+l i=l 

k m >0,Vm=l,...,d-h+l, 
k m =j~d+h 



Proof. 

1) The proof is equal to the proof of Lemma 3.2 in [2]. 

2) We can suppose e < d — 1, since otherwise the two members of 2) are equal to the 
null matrix. Set: 

r d-h+1 ^ 

= \ h,-,k d - h+ i : k m > 0,Vm = l,...,d- h+ 1, ^ £; m =j-d + /i>, 



/(A h ,...,Ad,i-d+/i) = J] n A 

-ffhj m=h 



m=l 

kd — m+l 

m 



Proceed by induction on j; if j = 0, V e = 0, 1, d — 1, the thesis is true. 
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Suppose that the thesis is true for j; then, for j + 1, Ve = 0, 1, d — 1: 



i=i V i=i / 

d d e 

= E { t A-\ h I + X h I)f(\ h ,...Ad,j-d + h) \{ (tA-X^UCA- Xil) 

h=e+l n=h+l i=l 

d e 

= ( t A-X e+1 I)f(X e+1 ,...,X d ,j-d + e + l) 11 CA-Xn^HCA-Xil) 

n=e+2 i=l 
d-1 / d 

+ E (*hf(X h ,...,Xd,3-d + h) Y\ ( l A-X n I) 

h=e+l V n=h+l 

d \ e 

+( t A-X h+1 I)f(X h+1 ,...,X d ,j + l-d + h) 11 ( l A - X n I) J JJ(*A - Xil) 



Observe that 

Xhf{Xh, -,Xd,j -d + h) + f(X h+1 , ...,Xd,j + l-d + h) = f{X h , ...,X d ,j + l-d + h), 

since f(Xh, X d ,j + 1 — d + h), a homogeneous polynomial of degree j + 1 — d + h in the 
variables X^, ...,X d , can be obtained by multiplying /(A^, X d , j — d+ h), a homogeneous 
polynomial of degree j — d+h in the variables Xh, X d , by the variable A/,, and by summing 
up /(A/j + i, X d ,j + l — d + h), a homogeneous polynomial of degree j + 1 — d + /i in the 
variables A/j+i, A^. Then: 



Y[C A - XJ) = /(A e+1 , Ad, j - d + e + 1) J](^ - XJ) 

i=i i=i 

cJ — 1 d e e 

+ £ /(A ft ,...,A d ,j + l-d + h) n (U-A n 7)n(^- Ai/J + A^JJCA- A,/) 

h=e+l n=h+l i=l i=l 

= e e n^- m+i n fA-A^ncA-v), 



n=/i+2 




+A^n( t A-A,7). 



1=1 



h=e+l Hh,j+l m=h 



n=h+l 



i=l 



since nf=i(*A — AjJ) is equal to the null matrix. □ 
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3 Main results 



Theorem 3.1. Assume that A has eigenvalues Ai, G C*, |Aj| ^ 1, for i = 1, 
and assume that the support of fi is not parallel to a proper subspace of Q fc invariant 
under A. Then, there exists c G R + such that, for all p G N such that gcd(det(A) , p) = 
gcd(det(-B),p) = 1, and for all n > c(lnp) 2 , we have: 

\\Pn — U\\ < e{p), where lim e(p) = 0. 

p— >oo 

Proof. From the lemmas 2.1 and 2.2: 

\\Pn-Uf<\ Y, l P «^)| 2 ' ^ 

aez£-{o} 

l^(«)| 2 = nV(^«)| 2 =nf E Mu)M(v)cos(-<u-v,^a>) J . (3) 

3=0 j=0 \u,vGZ fc V P ' j 

n—1 

In order to estimate J~J \\i ('A 7 a) | , V 5 € (0, |) set: 



3=0 



k 



a 



P 



L S = \j([o,i\ i - l x[s,i-8\x[o,i] k - i ), ^ = 
i=i 

and indicate by the vector whose components are the fractional parts of the compo- 
nents of £j. Consider the vectors yi, ...,yfc defined by Lemma 2.4; then, Vm = 1, k, 

y m = A Zm x m , where x m G V, z rn G {0, 1, d - 1}. 

Finally set z = max z m . The following results hold: 

m=l,...,k 

i+z 

Lemma 3.2. There exists b G (0, 1) such that, if {£J G L s , then J~[ |/2 (*A J a) | < b. 

3=0 

Proof. Consider the function g : [0, l] k — ► [0, 1] defined by: 

9 (t) = n I? c^pt) i 2 = n ( e ^( u )^( v ) c ° s ( u - v ' *^t» ) • 

j=0 j=0 \ u ,vgZ fc / 

If ti G (0,1) for some I G {1, ...,fe}, then pt ^ (mod p) and, by Lemma 2.4, 3 m G 
{1, k} such that (x m , l A Zm pt) / (mod p); then (x m , t A Zm t) £ Z and, by definition of 



g, we have g(t) < 1. Since g is continuous and L$ is closed and bounded, then g (t) < b < 1, 
V t G L5; in particular, if € 
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i+2 i— 1 2 

n |p (u'a) 1 2 = n 1? (*^'«) r n |# c^ +i «) " ^ ^ < &. □ 

i=0 j=0 Z=0 



Lemma 3.3. There exist <5,c G R + , (5 G (0, |), and j G N*, j < chip, such that, for 
any p G N sufficiently large and for any a G Zp — {0}, has a component in [5, 1 — 5]. 

Proof. The proof follows from Lemma 3 in [8] . □ 

Let c > — — % and suppose n > c(lnp) 2 ; then, for p sufficiently large, we have n > rt, 
mo 

where 

t = [clnpj + d, r = [clnpj , 

k 

for some c > — - — -. 

In 

Let j G N* defined by Lemma 3.3. If a G 7Jl — {0}, from Lemma 3.2 and Lemma 3.3, 



{%} G L ^ 



since 

j+z 



JJ |£ (U^a) f < 6 < 1. 
3=0 



Since j + z < t — 1, from (|3|): 



t-1 

|2 



iA(«)i 2 = ni^(^ a )r^ 6 - 

By repeating the previous arguments with t A tt a instead of a, V i = 0, 1, r — 1, we have: 

ip rt («)i 2 =nniMC^ +i *«)i 2 <^- 

i=0 j=0 

Then, from ([3]): 

E l^(«)! 2 < E \Prt(a)\ 2 <b r p k 
«ezfc-{o} oez*-{o} 

t£inp p k exp ((ilno+A;) hip) 



< 



b b 



ex P ((clnb+k) hip) _ 
Since lim = 0, by definition of c, from ([2]) we obtain the statement 

p — >oo b 

of the theorem. □ 



8 



Theorem 3.4. Assume that A has eigenvalues Ai, A& G C* such that A^* G N* — {1}, 
for some /j G N*, for any i = 1, and assume that the support of \x is not parallel to 
a proper subspace of Q fc invariant under A. Then, there exists c G R + such that, for all 
p G N such that gcd(det(^4),p) = gcd(det(B),p) = 1, and for all n > clnplnlnp, we have: 

\\Pn — U\\ < e(p), where lim e(p) = 0. 

Proof. Set 

I = lcm{k : i = 1, k), C = A 1 , crj = A-, V * = 1, k. (4) 

d d 
Since n - V) = 0, we have n (*C - = °- Moreover, VEC{1, d}, define: 
i=l i=i 

Yi? = | a G Z£ - {0} : ]J (*C-^/)a ^ (mod p), 
Yl( t C-a i l)a = (mod p), V T C {1, d} : \E\ + 1 < \T\ < d \ . 

i€T J 



The following relation holds: 



z£-{o}= (J Y * ( 5 ) 

In fact, observe that, V a G Zj - {0}, we can define: 

e = max i t G {0, 1, d - 1} : JJ (*C-<7jl)a ^ (mod p), 

for some -E $ {1, <i} such that = i} ; 

then a G Ye, for some E Q {1, d} such that = e, and this implies ([5]). 
Moreover: 

Y E = (J Y S)Ej 

0^5c{l,...,fc} 

where, if the vectors yi, ...,yfc are defined by Lemma 2.4: 

Y s ,e = | a G Y E : ^y m , JJ {'C-aj)^ + (mod p), V m G 5, 



y m , J] ('C-<7i/)a ) = (mod p), V m £ 5 1 . 
\ ieE / J 



Then: 



E ft( 



«;i 2 = E E Ei^(«)i 2 - 

a€Z£-{0} EC{i,...,d} 0^5c{l,...,fe} Vs.E 

If ^ S C {1, fc}, reorder the set S in the following manner: 

S = {mi,s, m|S|,s} , where m^g < mj,g i < j. 

Then, V/i = l,...,|5|: 

y™ h ,s = -4 Zm,i ' s x mhi5 , where x mft s G V, z m/l g € {0, 1, d - 1}. 



(6) 



Set x^g = x mh 5 , y ft|S = y mhjS , z KS = z mh S and let u^g, v hj g the vectors of the support 
of fi such that u^g — ^h,s = *h,S- From Lemma 2.2, V a £ Z* - {0}, we have: 



n-l 



\Pn(a)\ 2 < jj |£(*C"a) 



i=o 



^ II II ( 1 - Muh,sMvh,s) + 2/i(u fc) s)/i(v ft) s) cos f — (y h)S , *C J a) 



(7) 



where M h)S = { (h - 1) 



|5| 



n-l 



-^,5-U,V/ l = l,...,|5|. 



If a 6 for some % ^ S C {1, ...,&} and E 1 C {1,...,^}, reorder the numbers 

oi, ...,0<i so that the first | _E7 j correspond to the set {0^ : i <G E}. Moreover, V j, n <G N, 
V/i = set: 



\E\ 



ah,S,E = ( y h ,s, Hi'C-a^a \ E Zip — {0}, ^,n,h,S = ( 7 h ,s? C^C-aj) 



i=i 



i=i 



Q: 



P 



(the fractional part of ^y^g,* C J JJ(*C— 0jl) — y ). Observe that, from Lemma 2.5, 1): 



|E| \ |£|-1 



fj+|S|,0,h,S = \ { YKsf&Yl ( tC - <J i I ) ~ ) + E (7 m+lCj+|B|-m-l,m,/ l ,5 

i=l ^ / m=0 



> . 
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Moreover, use Lemma 2.5, 2) and the definition of Ye\ Vq £ Yg Ei V j G N, in the right 

.\E\ 

member of 2) multiplied by a (with e = \E\, A = C and Aj = cTj), only <x^n yC—(XiIja is 



different from (corresponding to h = m = d, k\ = j). Then, V h = 1, |S|: 

IE / 



\E\ \E\ I |S| \ 

i=l i=l \ i=l ^ / 

|E|-1 

j ah,s,E 



P 



j+IEI.O,/!,^ 



P 



—m—l,rn,h,S 



(8) 



m=0 



Let gh,s ■ [0, 1] — ► [0, 1] the function defined by: 

9h,s(t) = 1 - 2^(u h>s )lJ-(vh,s) + ty(uh,s)Kvh,s) cos(27rt). 
From ([7]), we have: 



\s\ 

^\p n (a)\ 2 < yi n n ^s^o, 

Y S ,E a htS . E eZ p -{0},h=lj€M hi s 
V h=l,...,\S\ 

\S\ 

= n n 9h,s (Cj,o,h,s) • 

fe=i afe,s,i5eZp-{o} jeM hiS 



h,s 



(9) 



Set L 



1 



1 



where cr = max <7j. 

i=l,...,fc 



2 d cr d+1 ' 1 2 af er< i + 1 

Observe that ^5 is continuous, gh,s(t) = 1 •<=>• i £ {0, 1}; then, 5 has a maximum 



6^5 < 1 in L, since L is closed and bounded; in particular: 

9h,S (£,j,o,h,s) < < 



&ft,s if £,j,o,h,s £ £ 
1 otherwise. 



(10) 



The following result follows: 



Lemma 3.5. Suppose £j,o,h,S,€j+i,o,h,S, -,Zj+e-l,o,h,S & L ' for some 3 £ N*, e £ N; 
then, for any s £ {0, 1, e — 1} and for any r G {j,j + 1, j + e — s — 1}: 

1 1 



€r,a,h,S £ 



2o!— Sg-d— s+1 ' 2°! — s (jd— s+1 



(mod Z) 



(11) 



In particular: 



sj+e-s-l,s,h,S £ 



1 



1 



2d-s a 3' 2 d - s cr3 



(mod Z). 



(12) 
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Proof. Prove the lemma by induction on s; by hypothesis, if s = and r £ {j,j + 
1, j + e — 1}, the thesis is true. 

Suppose that the thesis is true for s = n; then, for s = n+1, V r = j+e—n—2, 
we have: 

n+l n n 

i c r n (*c - = n (*c - a*!) - o n ^c- n cc - ^ . 
i=i 

By the inductive hypothesis: 



i=l 



i=l 



€r,n+l,h,S £ 2cr 



1 



1 



^d—riQ-d—n+l ' ri(jd— n+l 
1 



(mod Z) 



1 



2<i-(7i+l) ( j<i-(n+l)+l ' 2rf-(n+l) cr (i-(7i+l)+l 
Thus, we have (llip . In particular, since d — s > 2, (I12p follows. □ 



(mod Z). 



In order to finish the proof of Theorem 3.4, we will borrow some arguments from the 
papers [2], [3] and [6]. 

Fix h, S, E, set a = a^^E and consider the expansion of 2 in base cr^: 



P 



0. 010203... Define: 



* = \ l og ad P] (then, <r* 1 < p < o* d ), 



r = rh S 



\M h ,s\ 








- Zh,s\ j t 


t 




[( 


\s\ 



Moreover, recall that a "generalized alternation" between two consecutive digits OjOj+i 
of the expansion is defined as either the case Oj 7^ aj+i or the case a,j = a J+ i ^ {0, — 1}. 

Lemma 3.6. Suppose that there is a generalized alternation between the digits j + 
1, j + 2 - th of the expansion of |; then (,j+ifi t h,s £ for some i £ {0, 1, \E\}. 
Proof. The assumption imply: 



1 



a 



2 • 



1 



C 



-1 i_J_ 

a 2 ' a 2 
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If £j,o,h,S,€j+l,o,h,S, ->£j+\E\-l,o,h,S & L > then ' from © and Lemma 3.5: 



£,j+\E\,0,h,S G ( 


' 1 


1 - 


1 " 


+ 


< 


" 1 


1 - 


1 " 


+ 



\E\-\ \E\-\ 

E°Vn+l \ -* 

m=0 m=0 



m=0 



a* 2 d - m 

m=0 



1 



1 



1 



2 d CT 2 ' 2 d a 2 



c L. a 



Consider the first rt integer numbers of the set M^ g an d partition such numbers into 
r disjoint sets Mj = Mi^s, 1 < * < r, each of length t, such that, if % < j, x G Mj and 
y G Mj, then x < y. Moreover, V i = 1, consider the block of digits 

B a ,i = B a)i)h: s = {aj : j £ Mi} 

and, V B C {ai, 02, •••}, S made up of consecutive digits, indicate by A(B) the number of 
generalized alternations in B. Finally, VflcN, set: 

C{D) = C(D, h,S) = \{jeD: ^ fiAS G L}\. 

Suppose A({aj : j G M^s}) = m d, for some m G N. Since \E\ < d — 1, from Lemma 
3.6 we deduce C(Mf lt s) > m; in general: 



C{M KS ) > 



A{{ aj : j G M h , s }) 11 MB ' 



d 

It is possible to prove the following two results: 



> 1. (13) 

a 



Lemma 3.7. V i G {1, ■~,r}, as a ranges in Z p — {0}, the blocks B a ^ are distinct and 
have at least one generalized alternation. Moreover, V i, j G {l,...,r}: 

{B a>i : a G Z p - {0}} = {B aJ : a G Z p - {0}}. 

Lemma 3.8. 

3=1 »=1 J'=l 
13 



where, V i = 1, r, V j = 1, s, 7Tj is a permutation of {1, s} and a, > 0. 
By utilizing the formulas ([9]), (1101) . (I13j) and the lemmas 3.7 and 3.8, we have: 

Ei^Hi^ft E II 9 h ,s(ko, h ,s)<U E C /h ' s) 

Vs,s fi=l aSZ p -{0} j&M hiS h=l aeZ p -{0} 

, \S\ r \S\ 

< ^n e iut*" 

fc=l a£Z p -{0} i=l /i=l 



(where 65 = 6^5, /^g = y/b^ < 1) 



< 



, |S| |S| 

e /ir-"<^n e /; 



rA(B) 



5 11 ->n,S 

h=l length B=t 
A(B)>0 



h=l aeZ p -{0} 

Indicate with M(J) the number of blocks of length t with A(B) = j; then: 



(14) 



M(j) < 



t-i 

f rA(B) 

.,S 

length B=t j=l 
A(B)>0 



E C 



3=1 \ 

ffd ( (1 + cfdfh,s) 1 - l) < ( ex P {°dtf r h: s) - 1) • 



Suppose c/j 5 > 



/I5I 



3 Cfc « G 



w > s In pin In p. Since £ 



lnp 
In cr d 



(15) 



< ^ + 1, 



In cr d 



In f h S In CTrf 



,c/is such that, for p sufficiently large: 



\s\ 



ZhS / t 



hia d ch,s , , . . 2 hip 



trdtf^s ^j^- ex P ( ( 1 + TT^j ) lnln ^ ) = TM(P), 

where lim Jhsip) = 0> D Y definition of 5. Finally: 

p— >oo ' 



(US]) < cr d (exp (7 ft) s(p)) - 1) = e h S (p), where lim e ft) s(jp) = 0. 

p— »oo 

From the formulas ([2]), ([B]), (|14p and (|15|) . we obtain the thesis of the theorem, with 
c = maxc/,^. □ 
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Theorem 3.9. Assume that A has eigenvalues Ai, A^, A/% G C* such that A^ G 
N*, for some /j G N*, for any i = l,...,k, where |Ai| = 1 and |Aj| > 1 V i = 2,...,d, 
and assume that the support of fj, is not parallel to a proper subspace of Q fc invariant 
under A. Then, there exist a,c G R + and N G N such that, for all p G N such that 
p > N, gcd(det(^4),p) = gcd(det(B),p) = 1, and for all n > cp 2 , we have: 

|P n -C/||<2^expf-^ 



p 2 J 

d 

Proof. Define I, C and cij, V i = 1, k, as in (jlj); then, Y\ ( C ~~ ^i-O = 0. Suppose 

i=l 

d > 1 and define, V{l}cfiC {i ; d}: 

Z = |a G Zp — {0} : (*C-l)a ^ (mod p)} , 

Z E = la e Z :Yl (?C-<nl)a + (mod p), 

Yl (*C-CTj/)a = (modp), V T C : \E\ + 1 < |T| < d \ . 

i€T J 

The following relation is analogous to (|5]): 

Z= |J Zjj. (16) 

{l}CSC{l,...,d} 

Then: 

£ |P„(a)| 2 = £ ^|P n (a)| 2 + ^|P n (a)| 2 . (17) 

aeZ£-{0} {l}cBg{l,...,d} ^ ^ 

Analogously to the proof of Theorem 3.4, where we have valued |P n (a)| 2 , by utilizing 

Y S ,E 

the fact that > 1, V i = 2, d, we can prove that 3 c\, C2 G R + such that: 

^ ^ l-Pn(a)| 2 < £i(p,n), where ei(p,n) = ci In pexp ( ^— J . (18) 

Indeed, we must estimate only the sum |P n (a)| 2 . Observe that 

Z c = |J Z S , 
0^5c{l,...,fc} 

15 



where, if the vectors yi, - - - , are defined by Lemma 2.4: 

Z s = {a G Z c : {y m , a) ^ (mod p), V m G S, 
(y m , a) = (mod p), V m £ S} . 

Then: 



£|P n (a)| 2 = Yl El P «(«) 



(19) 

Z c 0^Sc{l,...,fc} Z s 

If ^ S C {1, A;} and /i G {1, IS 1 ]}, define x h>s , y h S , z h> s, Tih,S and v h ,s as in 
the proof of Theorem 3.4. Moreover, set: 

ah,s = (y~h,s, a ) ■ 
Observe that, V j G N and V a G Z c , we have: 

(y~h,S, tc3a ) = a h,S (mod p). 
Then, from the formula ([7|), V a G Z c : 

iPn(«)i 2 < n iMfc j «)r 

i=o 



< IJ ^1 - ^K u h,s)K v h,s) + 2 K u h,s)^h,s) COS ^ — a^s 



2^ 



(20) 



where |M^g| = 
Yl\Pn(a)\' 

Z S 



I n ~ 1 
L j 

|5| 



n — 1 , , . n — fciZfe + 1) 
2fc,5 > - 1 - (d - 1) > \. ; - Then: 



l\S\ 



Ik 



\s\ 



a hiS £Z p -{0},h=l 
V h=l,...,\S\ 

\S 



- II ( 1 ~ 2 K u h,s)K v h,s) + 2At(Ufe ) 5)/i(v fe) 5) cos 



2tt 



/ /2ir \\\ M h,s\ 

] J E ( 1 ~ 2 Ku-h,s)v(vh,s) + 2Ku-h,s)v(vh,s) cos ( — a hiS j j . (21) 

h=la h , s gZ p -{0} V 

2 2 

Note that —1 + cosx < ^x , V x G [— 7r, tt]. Furthermore, if a/^s G Z p — {0}, we can 



suppose: 



a h ,s € Z* n 



p — 1 p 
2~ ' 2 



2tt 



G [-vr,vr]. 
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Then, V h = 1,...,\S\: 

( /2-7T \\\ M h,s\ 

E ( 1 - 2 n(uh,s)Kvh,s) + 2Ku-h,s)v(vh,s) cos ( — a>h,S J J 
a h , s ez p -{o} V V P JJ 



< 2 E ( 1_ Zpt*( lt h,s)li(yh,s)<Xh,s) 

a h , s e(Zn[l,§]) V 

<2 E ex p(^-j^v(u h ,s)Kvh,s){n- (dlk + l))al jS J . 



r 3 1, 



Let r G (0, 1) such that 2r 3 + r 2 - 1 = (<£> — = and set: 

1 7" 2 

16 .„ V A -IkhiT 



th,S = exp (-jj^pfi(uh,s)Kvh,s){n - {dlk + 1))^ , c = max 



h,s 16n(u h>s )fj,(v h)i 
Let c > c, p sufficiently large and n > cp 2 ; then: 

n - (dlk + l)>cp 2 th,s < t, 

from which 

d22D = 2 ( e 4f J = 2 k*+ e *h$) z 2 E is s 

VlsSN* / \ a h ,s>2 / \ a h ,s>2 

= 2 ( tft ' s+ i~^) - 2 ^' S ( 1+ T^) = 3t/i ' 5 ' 

by definition of r. Then: 

/ /2vr \\l M fe,sl 

E ( 1 - ty(uh,s)l*(vh,s) + 2fJ>(uh,s)li(vh,s) cos ( — a fti s J J 
a, iiS ez p -{o} V 

< 3exp ^-|^/i(u ft) 5>(v fo)lS f)(n - (dZfe + l))j 



II E ( 1 ~ 2 m(u/ 1 ,s)m(v/ 1 ,s) + 2,u(u /ti s)/i(v/ li s) cos ( — a feiS J J 

ici / 2a s (n- (<«*; + , 8^1 , , , , 

< 3 1 1 exp 2 ' where a s = — > M u /i,sM v h,.s)- 

\ V J «* ft=1 
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Moreover, from (PT9l) and (|2T1): 



|S| 



(4^l) m p(- 2a '"-<f +1 >») 



where a = min = ^min//(uft j 5)/x(v/ lj 5) = — Then, from ([2]), (fT7|) and (fT8j) : 

0cSc{i,...,fc} /ijS 

II n rrll2 / 1 / , x /,* ,\ ( 2a(n - (dlk + I)) 
\\P n - U\\ 2 < - (ei(p,n)+ (4 fc - lj exp ^ ^ ~ 

2a(dlk + 1) 9 eAp,n) 

Observe that lim = and, if n > cp , hm — — = 0. Then, if p is 

p^oo p z p^ooexp (— 2an/p 2 ) 

sufficiently large: 

\\P n ~ U\\ < 7 4* exp 



from which 



p2 



\P„ I il <2 fc - 1 exp(-^] . □ 



Theorem 3.10. Suppose that the matrix A has an eigenvalue A G C, |A| > 1 (hence, 
the matrix l A too), that the support of \i is not parallel to a proper subspace of Q k 
invariant under A and that HB^H G L 2 , for all n G N. Then, there exist c € R + and 
N G N such that, for all p G N such that p > N, gcd(det(A),p) = gcd(det(l?),p) = 1, 
and for all n < chip, we have: 

\\P n — U\\ > —TjOp), where lim rj(p) = 1. 

2 p-^oo 

Consequently, O(lnp) steps are needed to reach the uniform distribution. 
Proof. If 1 A has an eigenvalue A G C, |A| > 1, then: 

ll*Ax|| 

11*41^= sup " "°° > |A| > 1. 

°° xGC ft -{0} ll x lloo 

Let a G Zp — {0}; from the lemmas 2.2 and 2.3: 

Pnipt) 



\Pn-U\\>\ 



n-1 



1/2 



ij] E Mh)Mi)cos(-(h-i,'^a) 
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X 2 

Since cos x > 1 — -ViGR,we have: 

1 n - 1 ( »||'^|| 2j \ 1/2 
\\Pn-U\\>-X\ 2 "°° , where = 2vr 2 A; 2 ||a||^ £ ^i(h)/x(i) ||h — G R" 1 

j=o V P / h,iez* 



(23) 



Moreover, 3 d G R + such that 1 — x > exp(— 2x), V x G [0, d]. 

Suppose n < chip, where c < — rn-rr, — . Then, V j = 0,1, ...,n — 1: 
" In Halloo 

l|U||i < HUII?. < ||U||!! nP =/HNIL < p 2(cln|N|oc- 1 



t 4 |l 2 ' 
F 



nun 2 "'' 

Since limp 2 ( cln ll ^Hoo -1 ) = 0, for sufficiently large p we can suppose G [0, d]; 

p^oo p 1 

hence: 

i n ~ l / „ IU A II 2 -? \ i / ™ _1 

ll-Pn-^H > o 1 1 eX P 







(-";"*) 


= r xp ( 



j=0 \ y J V 1 j=0 



1 / o W'Afl-Y 



■ exp , , 

2 V NiL-i p 2 J 

^ 1 I ^ 2(cln|| t A|| -l) , P \ _ 1 / n 

2 \ ll^llL-i (ll^llL-i)p 2 / 2 

By definition of c, we have the thesis. □ 

Theorem 3.11. Suppose that the matrix A has an eigenvalue A G C such that X 1 = 1, 
for some I G N* (hence, the matrix l A too), that the support of ji is not parallel to a 
proper subspace of Q k invariant under A and that HB^I^ G L 2 , for all n G N. Then, 
there exist 7 G R + and N G N such that, for all p G N such that p > N, gcd(det( J 4),p) = 
gcd(det(2?,p) = 1, we have: 

1 f l n 
2 eX H"F 



|Pn-^||>-exp(-^). 



Consequently, O (p 2 ) steps are needed to reach the uniform distribution. 

Proof. The assumption on A implies Ux = x, for some x G C k — {0}, and so 
(U - J)x = 0, which implies x G Q k - {0}; then 3 a G Z k - {0} such that l A l a = a. 
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y P > !l a lloo> we can suppose a G Zp — {0}, l A l a = a (mod p); then, V j € N, 3 i € 
{0, 1, I - 1} such that l Ma = t A i a (mod p). 

By proceeding as in the proof of the previous theorem, we obtain the following formula, 
analogous to (|23j) : 

i / \ n / 2 

\\P n -U\\ > - [ 1 



2 V P 2 , 

where 7 = 2-7T k I a „ max A > u(h)u(i) h — i| „ € R . Finally, V 

" lloo ie{o,i,.-,^i} " "°° . 



h,ieZ A ' 

sufficiently large: 

1 / "fn 
2^{-p 



P n -U\\>-exp(-lj).n 



4 Problems for further study 

A natural problem to study is the generalization of the recursion ([T]) to the analogous 
recursion in R fc reduced modulo c, for some real number c. In this case, the idea is to use 
the Fourier transform on R fc instead of in Z fc and then generalize the lemmas in Section 2. 
The expectation is to prove the convergence in law of the Markov chain and to estimate 
the rate of convergence to the uniform distribution on some subset of R fc : the set where 
the chain ranges. However, this set can be different from R fe (mod c) and it can be also 
countable. In order to establish it and to develop the theory, some changes of the results 
in Section 2 are needed: for example, Lemma 2.1 (upper bound lemma) is not valid in 
the continuous context and it must be modified. Conversely, results as Lemma 2.5 seem 
useful also in the study of the high powers of the real matrix A in the modified recursion 

(ID- 
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